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Abstract: We study the heavy charged Higgs boson (from 800 GeV to 1500 GeV in 
this study) in production associated with a top quark at the LHC with the collision energy 
yfs = 14 TeV. Such a heavy charged Higgs boson can dominantly decay into a top quark and 
a bottom quark due to its large Yukawa couplings, like in MSSM. To suppress background 
events and to confirm the signal, we reconstruct the mass bumps of the heavy charged Higgs 
boson and the associated top quark. For this purpose, we propose a hybrid-R reconstruction 
method which utilizes the top tagging technique, a jet substructure technique developed 
for highly boosted massive particles. By using the full hadronic mode of pp — > H^t — > ttb 
as a test field, we find that this method can greatly reduce the combinatorics in the full 
reconstruction and can successfully reduce background events down to a controlled level. 
The sensitivity of LHC to the heavy charged Higgs boson with two b taggings is studied 
and a 9.5a significance can be achieved when m^± = ITeV. 
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1. Introduction 

Probing the Higgs sector and unraveling the miracle of electroweak symmetry breaking 
(EWSB) is one of the most important goals at the Large Hadron Collider (LHC). In the 
Standard Model (SM), the Higgs sector includes a physical neutral Higgs boson and three 
would-be Goldstone bosons, while in most of extensions of the SM more than one physical 
Higgs bosons are predicted [1]. For example, two Higgs doublet model (THDM) extends 
the SM by adding one additional SU (2) doublet and results in a rich phenomena in Higgs 
sector by boasting two CP-even Higgs bosons, h and H, one CP-odd Higgs boson A, and a 
charged Higgs boson . The discovery of neutral Higgs bosons and the measurement to 
their properties could undoubtedly help people understand the EWSB mechanism, while 
the confirmation/exclusion of the existence of a charged Higgs boson at the LHC can 
pinpoint down the whole Higgs sector and further distinguish models. 

Experiments have invested quite an effort into the search of charged Higgs bosons. 
There is a direct limit of Mjj± > 78.6 GeV from the LEP searches by its exclusive decay of 

— > tv and — > cs [2]. At hadron colliders, the search approaches for a charged Higgs 
boson differ in term of its mass range. For the low mass range m H ± < nit, the signal for a 
charged Higgs boson is from top quark decay t — > H + b followed by the decay H + — > fu. At 
the Tevatron, due to the effects of phase space suppression for a heavy charged Higgs boson 
production, the search is mainly focused on the low mass range mjj± < uit- Currently there 
is no charged Higgs signal detected from top decay nor obvious deviation from Wtb vertex 
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measurement, which can put a constraint to THDM on the small and large tan /3 regions for 
a charged Higgs boson mass up to ~ 160 GeV [3]. However, due to its large collision energy 
of the LHC, the search for a heavy charged Higgs boson can be feasible via the gb — > tH~ 
production up to a large mass range [4] . For the large mass range mjj± > mt , the signal is 
from the main production process the gb fusion (gb — > tH~) followed by its decays. When 
a charged Higgs boson is heavy enough, its main decay mode can be H~ — > tb while other 
modes (say H~ — > tv) are small and can be neglected. 

When the heavy charged Higgs boson is around ITeV or higher, the top quark from 
its decay can be highly boosted and the detectors of LHC may not resolve all jets from 
its decay when the angle separation parameter is fixed to a specific value (say R = 0.5 
in anti-kt jet algorithm). Motivated by new technique of "jet substructure" [5, 6, 7, 8] 
developed for highly boosted massive particles, we propose a "hybrid-R reconstruction 
method" and investigate the full hadronical decay channel of the heavy charged Higgs 
production. Previous studies on this process focused on its semi-leptonic decay, i.e. tH^ — > 
tib — > jjblubb and mainly depend on triple-6-tagging [9, 10, 11] or four-6-tagging [12] to 
suppress the background. To our best knowledge, this is the first attempt to the full 
hadronically decaying channel for process pp — > tH^ — > tib. Similar to its semi-leptonic 
decay, the full reconstruction of its full hadronic decay mode may be difficult due to its 
large combinatorics and the large backgrounds from M+jets and QCD multi-jets processes 
[4]. Furthermore, the b jet from the boosted top may be difficult to be tagged. Therefore, 
we will use the top tagging and the b tagging for other isolated b jets as well as the full 
reconstructed objects in the final state to suppress the background. 

In order to identify objects in the signal, we reconstruct all objects (i.e. two W bosons, 
two top quarks, and a charged Higgs boson) in the signal. In our hybrid-R reconstruction 
method, we first use a larger cone size for the top tagger to identify a highly boosted top. 
Then we utilize the fact that the b jet from the charged Higgs boson can be reliably resolved 
by detectors with a smaller cone size and should be the most energetic jet to reduce the 
combinatorics in the full reconstruction. The remaining task is to identify the associated 
top quark which is always produced near threshold of the process where all jets from this 
top quark decay can be resolved by the smaller cone. In short, our hybrid-R reconstruction 
utilizes the advantages of both a larger cone (for a top tagging) and a smaller cone (for 
other resolvable jets). And it is found that this method can work quite well when the 
charged Higgs boson mass is above 800GeV. 

To be more precise to define our hybrid-R reconstruction method, we firstly adopt a 
large cone size in the top tagging method in order to capture the highly boosted top (John 
Hopkins top-tagger 1 [13] in Cambridge- Aachen algorithm with R optimized for each mjj± 
value). Then we adopt a smaller radius R (anti-fct with i?=0.4) to resolve jets in the event, 
which is efficient to reconstruct the associated top with a small transverse momentum and 
the isolated b jet from the heavy charged Higgs boson decay. We introduce a few crucial 
cuts to suppress background and adopt the neural network and boosted decision tree to 
optimize cuts. It is found that, with the help of John Hopkins top-tagger and hybrid-R 

It is a top tagging method proposed by John Hopkins group and we will introduce it in section 2.2 
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reconstruction method, the full-hadronic decay of a heavy charged Higgs boson at the LHC 
(y/s = 14TeV) is reachable. We also study the sensitivity of LHC to a heavy charged Higgs 
boson varying from 0.8TeV to 1.5TeV and find that LHC can cover this heavy mass region 
when the integrated luminosity is assumed to be 100 fb _1 . 

The remainder of this paper is organized as follows. In Sec. 2, we firstly give a brief 
review of jet algorithms and jet substructure. In Sec. 3, we introduce the "hybrid-R re- 
construction method" . And then we employ the "hybrid-R reconstruction method" armed 
with top tagging method to study the process pp — > tH^ — > tib with both tops in hadronic 
decay mode in Sec. 4. Our conclusions follow in Sec. 5. 

2. Brief Review of Jet Algorithms and Jet Substructure 
2.1 Jet Algorithms and Jet Substructure 

A Jet algorithm is to reconstruct a jet from the collimated spray of hadrons and to make 
comparison between data and theoretical prediction. Jet algorithms can be casted into two 
categories: the cone-based algorithms and the sequential recombination algorithms. The 
former one relies on the intuitive idea to put a cone along the dominant direction of the 
energy flow. The later one repeatedly recombines the closest pair of particles which look 
like undoing the QCD branching. In this paper, we will focus on sequential recombination 
algorithms including the kt [14, 15], the Cambridge- Aachen [16, 17] and the anti-fct [18]. 
For a recent review of jet algorithms and its application at hadron collider, we refer to Ref. 



Sequential recombination jet algorithms begin with a list of four-momenta of pseudo- 
jets which can be reconstructed from the track detectors and calorimeters, both electro- 
magnetic and hadronic ones, and recursively combine pairs of momenta into jets. The 
combination rules can be described as : 

1. Calculate the distance dij between all pairs of pseudo-jets and the beam distance diB 
for each pseudo-jet. 

2. Find the minimum in the list of dij and diB- 

3. If the smallest entry is dij, recombine i and j into a new particle and then return to 
step 1. 

4. Otherwise, if it is diB, define it as a final-state jet, and remove it from the original 
pseudo-jet list. Return to step 1. 

5. Iterate this process until the original list is empty, i.e., all pseudo-jets have been 
clustered to jets. 

The dij and diB in most popular sequential recombination algorithms for use at hadron 
collider can be parameterized as follows [18]: 



[5]. 



dij = mm(p t f,p t J) 
d iB = Ptf , 
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(2.1) 



(2.2) 
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where p is a parameter that is 1 for the kt algorithm, -1 for the anti-A^ algorithm and for 
the Cambridge- Aachen algorithm (CA algorithm). The kt algorithm clusters soft particles 
first, while the anti-fct algorithm recombines hard particles first. And the CA algorithm 
clusters by angle ranking. There are a few comments in order: 

• There is no standard jet algorithms. The jet algorithm and the parameters could be 
chosen flexibly. 

• All these algorithms are infrared-safe and collinear-safe. 

• The anti-fct algorithm gets circular jets area which act much like an idealized cone- 
based algorithm, while the kt algorithm and CA algorithm get irregular jets area. It 
is one of reasons why experimental groups favor the anti-kt from the viewpoint of 
detector geometry. 

• The kt and CA algorithms incline to combine the harder pseudo-jets or pseudo-jets 
with a large separation in the later, so it is expected that the combination in last step 
is the reverse of the main QCD branching or the 2-body decay of a massive particle 
which can be easily used in jet substructure study by iteratively tracing backwards the 
clustering history. However, the recombination history in anti-fct algorithm couldn't 
be utilized to related the branching and to get the subjets. 

The study of jet substructure has received considerable attention in recent years. At 
the LHC, the massive electroweak particles, i.e., W, Z, Higgs and top quarks, are always 
produced at a high transverse momentum (Pt 3> Mz)- In this case, their decay products 
tend to become collimated due to the large Lorentz boost factor. For highly boosted 
massive particles with hadronic decay, it is expected that a large cone size R could capture 
all the collimated decay products in a "fat jet". And the mass of this "fat jet" is one 
indicator of its origin [19, 20, 21]. However, high energy QCD jets can also generate large 
masses. In order to suppress QCD background to an acceptable level, one have to resolve 
more physical information in order to pick out the "fat jet" originated from the decay of 
highly boosted massive particles. One way is to anatomize the substructure and to identify 
the subjets in a fat jet to reconstruct the delicate kinematical variables. 

Along the line of this thought, the jet substructure technique has a boosting develop- 
ment in recent years including early studies in theory [23, 24, 25, 13, 26, 27, 28, 29, 30, 31] 
and in experiment [32, 33, 34, 35, 36]. For recent reviews on jet substructure, see Refs. 
[5, 6, 7, 8]. The earliest study of jet substructure is Ref. [23] in the context of a search for a 
heavy Higgs boson decaying to WW. More recently, the BDRS algorithm [25] developed a 
technique using the mass-drop and the filtering to transform the high-Pt WH, ZH(H — > bb) 
channel into one of the best channels for discovery of Standard Model Higgs with small 
mass at the LHC. The filtering [25], i.e., reclustering the fat jet of Higgs with refined small 
cone size Rfm, is supposed to further reduce the degradation of resolution on jets caused 
by underlying events. Besides filtering, some new techniques including "trimming", "prun- 
ing", "variable R" are also developed to refine the jet resolution [29, 30, 31, 37]. Application 
of the BDRS method to Higgs search in the MSSM context can be found at [38, 39]. 

Below, we focus on the top tagging algorithms. 
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2.2 Top Tagging 

Top quarks play an important role in understanding electroweak symmetry breaking and 
searching for new physics. Unlike the case at Tevatron where most of top quarks are 
produced near the threshold, at the LHC many boosted top quarks can be created. 

In the conventional method, top quarks are reconstructed by recombining three daugh- 
ters, i.e , three isolated jets in hadronic decay mode or one isolated lepton, one isolated jet 
and large missing energy in leptonic decay mode. And 6-tagging always is used to suppress 
the large QCD background. 

However, this conventional method confronts some problems for the top with high-i^. 

• The decay products become collimated and it is difficult to find three isolated daugh- 
ters to reconstruct the top when a cone size is fixed to a specific value like R = 0.5 
in the anti-Zc^ algorithm. 

• The 6-tagging efficiency drops significantly as the Pt increases because the tracks 
become too close to clearly identify the secondary vertex from the B decays in the 
silicon detector. 

• For its leptonic decay, it is a task to determine the missing energy in the hadron 
collision environment. 

So, the jet substructure technique used to identify boosted tops from its hadronic decay 
have been developed in recent years. As summarized in Ref. [5], many discriminating 
variables are used on subjets for top tagging, such as d^-types variables [34, 35, 40], -Pi- 
ratio type variables [13, 26, 28, 30, 31, 40], event-shape variables [41, 42], and constraints 
on a W-subjet mass[13, 26] as well as other subjets correlation variables ( a helicity angle 
9t [13] ). In addition, the HEPTopTagger[28, 43] is designed for moderately boosted tops 
and some techniques are developed for leptonic top tagger [44]. 

Here, we adopt the top tagging method for the highly boosted top in hadronic decay 
proposed in Ref. [13]. We will dub it as "JH top-tagger" for simplicity and use it in later 
analysis. As well-known, a boosted top from its hadronic decay looks like a fat jet with 
three hard cores. Similar to the BDRS jet substructure method for boosted Higgs, the JH 
top-tagger firstly uses a large cone to cluster the event in order to capture all the decay 
products and relevant radiations of a boosted top and then declusters the top jet to find 
three subjets. 

To resolve a fat top jet into the relevant hard substructure, the following recursive 
procedure are applied in JH top-tagger [13]: 

1. Four momenta of particles is clustered with a large cone size R ( CA algorithm are 
used in original paper) 

2. The last cluster is undone to get two objects j% and j%. If the Pt ratio of the softer 
jet j2 over the original jet j is too small, i.e., P^ /P tj < 5 P , throw the softer j2 and 
go on to decluster on the harder one. 
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3. The declustering step is repeated until two separated hard objects are found. If any 
criterion below are satisfied, the event is failed: 1) both objects are softer than S p 
(2) the two objects are too close, Ary + A<fi < 5 r (3)the original jet is considered 
irreducible. 

4. Declustering repeatedly on these two subjets will result in 2,3, or 4 hard objects. 

5. The case with 3 or 4 subjets are kept. And then require the these subjets that the 
total mass should be near mt, the mass of two subjets should be in the mw window, 
W helicity angle 6t should be consistent with a top decay 2 . 

The parameters involved in the method can be optimized event by event [13]. In our study, 
these parameters are fixed as below: 

5 P = 0.19 5 r = 0.1 \6 t \ < 0.65 . (2.3) 

An interesting application of the Higgs tagging method as well as the top quark tagging 
method for the heavy top partner search can be found at [45] . 

3. Heavy Charged Higgs in Full Hadronic Decay Mode 
3.1 Heavy Charged Higgs 

For a heavy charged Higgs boson (m#± > mt), the main production process is the gb fusion 
process gb — > tH~ . And the LHC can probe the charged Higgs boson in a wide mass range 
up to TeV scale. Additional productions including W ± H T [47], $i? ± ($ = A, h, H) [48] 
and H + H~ [49] have also been studied but the cross sections are generally small in these 
modes. For the with m H ± > rrit + m^, the decay H~ — > ib is the main decay mode due 
to the large Yukawa coupling. Although the decay mode H~ — > tv is an interesting one, its 
small branching ratio makes it difficult to probe a heavy charged Higgs boson. The early 
studies on the search for a charged Higgs on pp — > tH^ production can be found in Ref 
[9, 10, 11, 12]. The search in this channel are also carried out by both ATLAS and CMS 
collaborations. These collaborations mainly investigate gb — > tH~ — > tib — > bqqT{had)vbb 
and the semi-leptonical channel gb — > tH~ — > tib — > bqqblvb. For search in r channel 
gb — > tH~ — > tib — > bqqT(had)ub, the branching ration of — > tv is small especially 
for higher mass. The main backgrounds to this signal are ti events, VF+jets and 
QCD multi-jet events [4]. For semi-leptonically decaying channel of production, i.e., 
ti7 ± ->• tib jjblvbb, the main backgrounds come from tt + jets, tt + bb(QCD), and 
ti + bb(EW) [4]. These studies indicate the search for a heavy charged Higgs boson from 
this channel has to rely on triple-6-tagging [9, 10, 11] or four-6-tagging [12] in order to 
suppress the background. 

2 The helicity angle 8 t is defined as the angle, measured in the rest frame of the reconstructed W, between 
the reconstructed top's flight direction and one of the W decay products. The lower Pt subjet in the lab 
frame are chosen to set the angle. 
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Armed with the jet substructure method to tag highly boosted top quarks, we inves- 
tigate the full hadronic decay mode for process pp — > gb ->• t 2 H~ ->• t 2 tib^. To our best 
knowledge, this is the first attempt for such a study in literature. For convenience, we 
here call the top quark and the b quark from the charged Higgs boson decay as t\ and 
63, respectively, and the associated top quark as t 2 as well. We attempt to reconstruct 
all objects in the final states. For this purpose, we resort to the top tagging method and 
a "hybrid-R reconstruction method" which can help to overcome the combinatorics prob- 
lem in the full reconstruction while maintaining a remarkable suppression to backgrounds 
(especially QCD background). 

Since we only concentrate on the study of the interaction of the heavy charged Higgs 
with top and bottom quarks, we simplify the Lagrangian as follows: 

C = C SM + dnH^dW* + mjj^H^ 

+H ± t{Y L P L + Y R P R )b + h.c. . (3.1) 

The couplings Yl and Yr are treated as free parameters in the model-independent analysis. 
In two Higgs doublet model, for simplicity, we assume the branching ratio of charged Higgs 
into tb is 100% and take Y L = Y R = 1. 

3.2 Hybrid-R Reconstruction Method for Full Hadronic Decay Mode 

For the full hadronically decaying of process pp — > gb — > t 2 H~ — > t 2 t\b^^ the main back- 
grounds are tt and QCD multi-jets events. We focus on a heavy charged Higgs boson with 
mass around 1 TeV or so. 

In such a mass region, the daughter t\ from the heavy charged Higgs boson decay is 
always highly boosted. And a larger cone size R is expected to capture all the hadronic 
decay products of t\ in a fat jet. However, the associated t 2 in the pp — > t 2 H^- process is 
produced near its threshold and its decay products generally can not form a fat jet. So the 
large cone size used for t\ tagging fails to catch the t 2 and only a small cone size in a jet 
algorithm can resolve all its three jets. Then the traditional top reconstruction method by 
recombining three jets is expected to be efficient to tag t 2 . 

In order to reconstruct all objects in the process pp — > gb — >■ t 2 H~ — > t 2 iibs and to 
exploit the advantages of a larger and a smaller cone size jet algorithms, we propose a 
"hybrid-R reconstruction method" . 

The five steps we used in our "hybrid-R reconstruction method" can be described as 
follows: 

1. At the preselection level, we adopt the CA jet algorithm with a larger cone size 
parameter and use a top tagger method to capture a highly boosted top and preselect 
events. The JH top-tagger is employed in this work. The jets at this step form the 
jet set Jo- 

2. After capturing a top jet labeled as ti, we re-cluster the pseudo-jets in the event with 
a smaller cone size parameter (the anti-fct algorithm with R = 0.4 is chosen as an 
example) and then get a new jet list Lq. 

3 We also include pp —¥ gb —¥ i 2 H + — y t2tit>3 in our calculation 
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3. If a small-size jet is within the larger cone of the direction of the tagged top jet, we 
remove it from the jets list for the 63 and t<i reconstruction. The rest of jets form a 
new jet list L\. 

4. We identify the most energetic jet in the list L\ as a 63 jet. This choice also implies 
that the tagged top t\ and the identified 63 fly back-to-back. This is due to the 
observation that the tagged top in the larger cone size by the top-tagger method 
contains subjets which are massive and energetic and the leading two objects with 
large transverse momentum in the smaller cone size always fly back to back. So we 
don't impose a cut on A(f> or AR here. The unused jets form the list L<i. 

5. We require there are at lest 3 jets in the jet list Li in an event. Then we introduce a 

I f ' ' \ PDG 1 2 

X 2 method to tag the second top quark in the signal. Here, x 2 = Jl 2 mw — — l~ 

PDG 2 ° W 

j2j3^ m t — [_ w j iere aw anc f at are gq U£ j t 10 GeV and 15 GeV, respectively, 

which is to account for the detector resolution capability. The and to^ dg are 

taken as 80.5 GeV and 173 GeV, respectively. The combination which minimize the 
X 2 is chosen as the right reconstruction. 

In step 2 we recluster a selected event with a smaller cone size, and such a reclustering 
has also been used to reduce the degradation of the resolution on a fat jet from highly 
boosted massive particle in filtering [25] and some other grooming techniques, and also 
can be used as preselection rule [28], and can be used to separate top jets from QCD jets 
in the hollow cone method [46] as well. But it should be stressed that, in our hybrid- R 
reconstruction method, the reclustering and removing the smaller jets within the radius of 
the tagged top is to catch the highly boosted top and to reconstruct successfully the other 
massive particle with a smaller Pi (Here it is the t<i). 

Another remarkable point is that any top tagging algorithm for highly boosted top can 
be used in the first step although the JH top-tagger is used as an example. 

Before carrying on the detailed analysis, we would like to show some features of the 
signal which support our hybrid-R reconstruction method. At parton level, in Fig. 1, we 
show the distribution of transverse momenta of t2, t±, and 63, and the angle separation 
among two of t\, ti and 63 in the case of M H ± = 1 TeV. The distributions of -Pt(ii) 
and Pi (63) are similar and yield an important clue for Pt{t\) and P t {b^). We find that 
Pt{ti) = Pt(bs) > jQ r mH ± can work and so we will use it in our simple cut method which 
will be introduced in the following section. We also show the distribution of Pt(t±) — Ptit^) 
and R{t\b^) — R(t2b^), which provide important hints to determine the right combination. 
As shown in Fig. 1, one can find the fact that in most of phase space the transverse 
momentum of t± is larger than ti. Meanwhile, the angle separation between t\ and 63 
demonstrates that these two objects indeed fly back-to-back. And in most of cases, the 
angle separation between t\ and 63 is larger than that between and 63. At parton level, 
it is also found that the A.<fi(t\,bz) is always close to tt, which may be useful in the full 
reconstruction. 

At the fast detector simulation level by using PGS4 [50], as demonstrated in the left 
panel of Fig. 2, we show the mass distribution of the leading massive jet, which is assumed 
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Figure 1: The transverse momenta and the angle separation are exposed at parton level with 
m H ± — l.OTeV. The transverse momentum difference between t\ and t 2 and the angle separation 
difference between R(tib^) and ^(£2^3) are also shown. 



to be a t\ jet in our hybrid-R reconstruction method. Then from the rest of jets, we pick 
a jet with the largest transverse momentum as the 63 jet. After picking out these two jets, 
we examine the invariant mass of the combined Lorentz vector from them, which is also 
displayed in the right side of Fig. 2. We examine how these two types of mass distribution 
vary to the change of angle separation parameter R in the anti-A^ algorithm. 

In Fig. 2, we can read out that when R = 0.4, neither a W boson nor a top quark can 
be captured in a single massive jet significantly. When R = 0.7, a considerable fraction of 
W boson can be captured in one single massive jet while a small fraction of t quark can be 
captured. When R = 1.0 and R = 1.3, quite a fraction of top quark can be captured by a 
single massive jet. 

Meanwhile, when R is large enough, the simple reconstruction method by using the 
most-massive jet and the most-energetic jet can even reconstruct the heavy charged Higgs 
boson mass bump, as demonstrated in Fig. 2. Both parton level and fast detector simu- 
lation level analysis demonstrate that our event reconstruction method captures the most 
important feature of the signal. 

In our Monte Carlo study, the signal events are generated by Madgraph/MadEvent 
[51] and background events are produced by Alpgen [52]. These events are fed to DECAY 
to generate full hadronic decay final states and pass to PYTHIA [53] to simulate showering, 
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fragmentation, initial state radiation, final state radiation, and multi-interaction as well. 
After that, we use fastjet [54] and SpartyJet[55] to perform jet clustering and analysis. 




Figure 2: The mass distribution of the leading massive jet and the reconstructed charged Higgs 
boson mass are shown with the case mjj± = ITeV. The dependence on the angle separation 
parameter R in the anti-fc t algorithm is also exposed. To be more realistic, the detector effects have 
been taken into account by using PGS4. 



Now, we examine the top tagging efficiency of JH top-tagger for different values of the 
cone size parameter R in both kt and CA algorithms by using the SpartyJet package. The 
results are tabulated in Table f and Table 2, which are two cases study for m H ± = l.OTeV 
and rriff± = 1.5TeV, respectively. As shown in Table 1 and Table 2, the tagging efficiency 
of the CA algorithm is better than that of the kt algorithm in a larger R, especially for a 
larger value of m#± . Another remarkable point is that the efficiency of the CA algorithm 
does not depend upon the value of R as that of the kt algorithm does in large R range. 
Due to these feature, in later analysis, we will use the top-tagger in CA algorithm with 
optimized cone sizes for different mjj±. 
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Table 1: The efficiency of top tagger in k t and CA algorithms along with the variation of angle 
separation parameter R are tabulated. The charged Higgs boson mass is assumed to be 1.0 TeV. 

This behavior can be understood by the clustering sequence of jet algorithms. In the 
CA algorithm, the clustering is performed by angle ranking and the softer drop method, 
i.e., 6 P cut in the top-tagger will efficiently remove the softer branching, underlying events 
and pile-up with a large separation to the hard subjets. Therefore it is successful to find the 
refined radius for each hard core in a top jet. However, in the kt algorithm, the clustering 
recombines the hardest and largest separated particles in the last stage of recombination, 
the softer pseudojets have been merged ahead and the 5 P cut actually don't work as it 
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CA algorithm 


23% 
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30% 
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29% 



Table 2: The efficiency of top tagger in k t and CA algorithms along with the variation of angle 
separation parameter R are tabulated. The charged Higgs boson mass is assumed to be 1.5 TeV. 

does in the CA algorithm. Additionally, in the kt algorithm, a jet has a larger area and a 
stronger Pt dependence than that in the CA algorithm [56]. 

These features also explain why the efficiency of top-tagger in the CA algorithm with a 
smaller R is smaller than that in the kt algorithm. In a smaller cone size, the CA algorithm 
can not capture all jets from the top quark decay while the kt algorithm is better since it 
can bring more protojets and get a larger mass for its larger jet area and no softer drop. 

It is also remarkable that the configuration of the b quark and the W boson from 
t\ decay flying near to the direction perpendicular to the boost direction of t\ has larger 
probability to form a top jet. While the configuration either the b or the W boson flying 
parallel to the boost direction can not form a top jet. The parameters of top tagging given 
in Eq. (2.3) also affect the tagging efficiency. This explains why there are more than 70% 
or so in the signal can not be successfully tagged as a top jet. 

4. Analysis Results 

In this section, we present our detailed analysis. One obvious question is whether signals 
can be recorded at detectors of LHC. Below we use the triggering system of ATLAS detector 
as an example to examine whether the signals can be recorded. 

At the level-1 triggering, signals can be selected out by the jet triggering from heavy 
charged Higgs boson decay which have large Pt- The mutli-jet trigger should also help 
since there are more than 3 jets in the signals. Furthermore, the scalar sum of transverse 
energy, which should be close to the heavy charged Higgs boson mass, can also help to pick 
out the signals. At the level-2 triggering, the tagged b jets can be helpful to record signals. 
At the event filter level, the top tagger algorithm can be implemented to select signals. 

At the offline analysis level, we introduce a set of cuts and compare their efficiencies 
between the signal and backgrounds. Multivariate analysis methods, the neural network 
and the boosted decision tree, are also used to improve the significance. 

4.1 Signal and Backgrounds 

The cross section of the process pp — > gb — ^ tiH~ — > £2^163 and the decay width of 
heavy charged Higgs boson are shown in Fig. 3. The cross section is about 0.2pb for 
m H ± = 0.8TeV, which drops down to 0.02pb for mjj± = 1.5TeV. While the decay width 
is about 75GeV for m^± = 0.8TeV, which increases to 170GeV for m H ± = 1.5TeV. It 
is found that the decay width can be approximated by ^m H ± ~ 10% x m h ±. After using 
the PGS4 to simulate the detector effect, it is also found that the effect of decay width to 
the mass bump of is roughly comparable to the detector effect, which means the mass 
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Figure 3: (a) The cross section as a function of M H ± at the LHC with yfs = 14 TeV is plotted, 
(b) The width as a function of M H ± is shown. 
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Figure 4: For the case m H ± = ITeV, the normalized distributions of signal and backgrounds for 
kinematic observables, Pt(ti), P t (bs), and H t , are demonstrated, respectively. 



bump of the charged Higgs boson is not too broad to be detectable from the detectors of 
the LHC. 

Since the cross section of SM background processes is huge, it is necessary to use the 
6-tagging (assuming 6-tagging efficiency 60% and mis-tagging efficiency for g and quarks 
1%) to suppress the SM background events. We demand that the most energetic isolated 
jet and one extra b jet in the jet list Lq must be tagged as b jets. The most energetic 
isolated jet must be a tagged b jet can suppress the ti + jets background by a factor of 100 
since the isolated jet in this process is mainly from gluon radiation. 

Below we introduce a few crucial kinematic observables. One of the most characteristic 
features of the signal is that the t\ and 63 from the decay of heavy charged Higgs boson 
are very energetic and fly in the small r\ region. Therefore, the top-tagger is efficient to 
capture the highly boosted t\ coming from the decay of the heavy charged Higgs since the 
large mass of m H ± will transform to the high Pt of the top quark t\. For tij background, 
there is only a fraction of events could provide a high Pt top jet passing the JH top-tagger. 



- 12 - 
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Figure 5: For the case m H ± = ITeV, the normalized distributions of both signal and background 
for reconstructed mass peaks, mw 2 , mt 2 , and rn H ± are shown. 



However, it is difficult to find a massive jet has the same substructure as the highly boosted 
top in QCD multi-jet events. Meanwhile, the Ht observable, i.e, the scalar Pf sum of final 
states, should be a good discriminant for signal and background separation. 

After tagging a highly boosted t\, for the signal the 63 jet generally should be the 
most energetic jet in the rest of jet list L\. It is also expected that the distribution of 
transverse momentum of t\ and 63 should be good for signal and background separation. 
Once both t\ and 63 are identified, the mass bump of the heavy charged Higgs boson can 
be reconstructed. Then the only combinatorics is from the reconstruction of the second 
top quark, t2 , which can be reconstructed from the remain jet list L2. While it is difficult 
for backgrounds to reconstruct such a top quark especially for QCD multi-jet. So it is 
obvious that the reconstruction mass bumps from ti including mw 2 and m; 2 are important 
to separate the signal and backgrounds. 

In Figs. 4-5, distributions of kinematic variables of events which pass the preselection 
rules are presented. As shown in Fig. 4, the distribution of Pt(t\), Pt(bs), and Ht for both 
signal and background events are shown. In Fig. 5, the mass bump of mjy 2 , m t2 , and m H ± 
are shown. 

According to these kinematic features, a set of preselection cuts are chosen and the 
comparison of the cuts efficiencies between the signal and backgrounds are shown in Table 
3 and Table 4. These two tables correspond to the case for heavy charged Higgs with mass 
1 TeV and 1.5 TeV at the LHC ( yfs=14: TeV), respectively. Due to the preselection rules 
introduced in our reconstruction method, the more relevant QCD multijet background are 
generated by requiring the leading four jets with Pt^min 

> lOOGeV. 

Some of the preselection rules, which have been encoded in our reconstruction method, 
are listed as 

• Detector acceptance cuts are chosen as Pt(j) > 20 GeV, \r/\(j) < 2.5. 

• A cut on the number of jets in the jet list Lq, nj > 6. 

• A cut on the Ht, Ht > tqM h ±, and the centrality, C > 0.3. 
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• A cut on the leading energetic jet in the jet list Lq, E(ji) > ^M H ±. 

• A cut on the second leading energetic jet in the jet list Lq, E(j\) > \M H ±. 

At the data analysis level, for the sake of comparison with the multivariate analysis 
methods, we introduce a simple cut method consisting of several important cuts after the 
full reconstruction, which are listed below 

• A large Pt requirement on t\ jet captured by the top-tagger and on the supposed 63 
in hybrid-R reconstruction method, i.e. Pt{bs) > j^M H ± and Pt{t\) > ^qM h ±. 

• Cuts are chosen for the mass bumps in the t 2 reconstruction, i.e., \my/ 2 — m^ >G \ < 20 
GeV and \m t , 2 - m t PDG | < 30GeV. 

• The charged Higgs boson mass window is chosen as \m H ± — m^™ med | < 200 GeV. 

The efficiency of each cut is also shown in Tables 3 and 4. The effects of b taggings, which 
can help to reduce combinatorics and suppress QCD background events, are also shown. 
Here the b tagging requires the leading P t jet must be tagged, where we assume 6-tagging 
efficiency as 0.6 and the mistagging efficiency as 0.01. 





signal 
m H ± = 1 TeV 


tt + jets 


QCD 

n, >4 


Cross Section x Br (pb) 


0.053 


553 


9186.8 


After JH tagger & H T > 400GeV 


29% 


7% 


4.8 x 10" 3 


Two 6-taggings 


10% 


4.2 x 10~ 4 


4.8 x 10~ 7 


The number of jets in jet list Lq, rij > 6 


10% 


3.7 x 10" 4 


3.4 x 10" 7 


H T > 700GeV & C > 0.3 


9.6% 


2.1 x 10" 4 


3.0 x 10~ 7 


The leading jet E(ji) > 350GeV 


8.9% 


1.6 x 10" 4 


2.5 x 10~ 7 


The second leading jet E(j 2 ) > 250GeV 


7.9% 


1.3 x 10" 4 


2.2 x 10~ 7 


P t (b 3 ) > 300 GeV & P^h) > 300 GeV 


5.3% 


4.5 x 10" 5 


7.7 x 10~ 8 


\m W2 -mP? G | < 20 GeV 
& \m t2 -mP DG | < 30 GeV 


3.1% 


2.6 x 10~ 5 


2.3 x 10~ 8 


\m H ± - m|f ± umed | < 200 GeV 


2.5% 


1.5 x 10" 5 


1.5 x 10" 8 


Events in 100 ftr 1 


133 


830 


13.8 



Table 3: The efficiencies of both signal and background before and after using JH top tagger are 
displayed. The collision energy is assumed to be 14 TeV. The JH top tagger and CA algorithm are 
used with the angle separation parameter R = 1.1. 

As shown in Table 3, the top tagging efficiency is 30% for signal while the mistagging 
efficiency for QCD background are at 0.1% level. For ttj background, there is a small 
fraction about 7% where one of the top quark can be boosted. As for M H ± = 1.5 TeV, 
the top tagging efficiency and the mistagging efficiency are similar. The comparison of the 
cuts efficiencies for the signal and backgrounds are shown step by step. After imposing 
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signal 
mu± = 1.5 TeV 


tt + jets 


QCD 

rij > 4 


Cross Section x Br (pb) 


0.007 


553 


9186.8 


After JH tagger & H T > 400GeV 


29% 


3% 


3.3 x 10- 3 


Two fc-taggings 


10% 


1.8 x 10" 4 


3.3 x 10- 7 


The number of jets in jet list Lq, rij > 6 


9.3% 


1.7 x 10- 4 


2.5 x 10" 7 


Ht > 1050GeV& C > 3 


8.3% 


5.5 x 10" 5 


8.9 x 10 -8 


The leading jet E{ji) > 525GeV 


7.7% 


4.4 x 10" 5 


7.3 x 10~ 8 


The second leading jet E{j2) > 375GeV 


6.8% 


3.4 x 10~ 5 


5.4 x 10~ 8 


Ft (63) > 450 GeV k P t {ti) > 450 GeV 


4.4% 


9.2 x 10" 6 


8.5 x 10~ 9 


|mw 2 -m™ G | < 20 GeV 
|mt 2 -mf DG | < 30 GeV 


2.5% 


4.6 x 10" 6 


4.3 x 10~ 9 


\m H ± - m^j umed | < 200 GeV 


1.8% 


1.9 x 10~ 6 


1.8 x lO" 10 


Events in 100 ftr 1 


13 


105 


0.2 



Table 4: The number of events of both signal and background before and after using JH top tagger 
are displayed. The collision energy is assumed to be 14 TeV. The JH top tagger and CA algorithm 
are used with the angle separation parameter R = 0.9. 



the requirement of successful reconstruction of £2 > the background have been considerably 
rejected. Furthermore, if a mass window cut on the reconstructed charged Higgs boson 
is imposed, the background is significantly suppressed by a factor 3 or so for the SM 
background while the signal is retained in 50% (for M H ± = ITeV). It is of interest to 
note that the cuts used in above steps can suppress the background by a factor 10 (100) 
for m.jy± = ITeV (mjj± = 1.5TeV) at least. This suppression factor increases with the 
increase of charged Higgs boson mass. 

We would like to make a comment from Table 3 and Table 4. These two Tables provide 
crucial evidences to use the top tagging substructure technique and b taggings. Without 
the top tagger and b taggings to select the most relevant events, the computing and analysis 
time for background events will increase drastically, which might be difficult to make the full 
hadronic decay mode realistic. Furthermore, if the standard top reconstruction methods 
are used for t\ reconstruction, we have to confront with the combinatorics problem. Even 
the most energetic top quark t± can be successfully reconstructed, it is found that the 
dominant background ti+jets can not be sufficiently suppressed. Similarly, without the 
top tagging, the mistagging rate from QCD background can reach to a few percent. For 
the sake of comparison, we have considered a full reconstruction without using the top 
tagger. We have considered events with number of jets rij > 7 and have tried a biased x 2 
method, where x 2 1S defined as 



2 _ \m Wl - m ™Z\2 | mw2 _ m Pp Q |2 Ki _ m PDG|2 | mt2 _ m P DG|2 
A — „2 T o "I" o, -T „ 2 



'W "W "t "t 

-|2 
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In order to reduce combinatorics and to save processing time, among the leading two jets 
of each event, we assume one must be 63 and the other must from t\ decay. The last 
term in Eq. (4.1) is to capture the configuration with t\ and 63 from the heavy charged 
Higgs boson decay, and we have used = 0.2. Then the mass of the heavy charged 
Higgs boson can be reconstructed. Without the help of TMVA methods, for m H ± = 1 
TeV (mjj± = 1.5 TeV) case, we find the significance can not be better than 4 (1.5) due 
to a large number of background events and the deterioration of wrong reconstruction of 
signal events. While when top tagger is used, we can greatly reduce combinatorics and 
analysis time and maintaining a better significance. Therefore the top tagging methods 
adopted here are indeed necessary. It is also remarkable that the b taggings are also crucial 
to suppress the SM background (especially the QCD background) to a computable and 
controllable level. Without b taggings, the significance can not be better than 1.2(0.5). 




Figure 6: The lego plot of Pt and 2-Dimension plot of 77 — $ for a typical signal event clustering 
with anti-fct R = 0.4 for illustration. And the direction of the top jet tagged by JH top-tagger and 
the top area with R = 1.0 are marked with a red circle. 



In Fig. 6, we also show the Pt lego plot and 2-Dimension rj — (ft diagram of a typical 
signal event (anti-fct R=0A) passing our selection rules. The direction of the top jet tagged 
by JH top-tagger and the top area with R = 1.0 are marked with a red circle. In this 



example, the traditional recombination method fails because only two jets can be found 
with the jet algorithm parameter R=0A within the tagged top jet. However, such a signal 
event can be efficiently selected out by our hybrid-R reconstruction method equipped with 
the top tagging algorithm. Furthermore, our method is efficient to reconstruct the other 
top because we reduce substantially the combinatorics after subtracting the jets within the 
tagged top and removing the 63 jet. For a naive comparison, the combinatorial number in 
the traditional recombination method is C| = 35 to reconstruct a top in one event with 
seven jets. In our method, the combinatorial number is Cf = 4. 

4.2 Multivariate Analysis and Search Bounds 

In this subsection, we explore the multivariate analysis (TMVA) to optimize the cuts and 
improve the sensitivity. Both of the neural network discriminant analysis and the boosted 
decision tree discriminant analysis are considered. TMVA methods have been broadly used 
in the top quark precision measurements by experimental collaborations [57, 58, 59] and 
W-jet substructure analysis [60], as well as new physics search studies, like Higgs boson 
search [61] and heavy quarks search [62, 63, 64]. The boost decision tree method has been 
adopted in the light charged Higgs boson search [65]. 




GeV GeV GeV 

Figure 7: For the case m H ± = ITeV, the normalized distributions of both signal and background 
for a few kinematic obscrvablcs, jj , and 77167 are shown. 

In Fig. (7), three extra kinematic observables are shown. Except the kinematic vari- 
ables introduced in Section 4.1, we also find several extra kinematic observables helpful to 
separate signal and background. 

• The observable rrij lt j 2 is the maximal invariant mass of a pair of jets in the jet list 

Lq. 

• The observable is the invariant mass of the leading two jets in the jet list Jo. 
This observable can form a bump for signals near the mass of heavy charged Higgs 
boson. 

• The observable 772-67 is the minimal invariant mass of a pair of jets for the leading 7 
jets in the jet list Lq. The signals tend to condense near the mass of W boson. 



-17- 



• The maximal value of invariant mass of leading three jets in the jet list Lq, which is 
labeled as mf 2 3) is found to be concentrated near the mass of heavy charged Higgs 
boson for signal. 

• The minimal value of invariant mass of three jets among the leading 7 jets in the jet 
list Lq, which is labeled as m^j, is found to be concentrated near the mass of top 
quark for signal. 

• The invariant mass of the leading 7 jets in the jet list Lq is found to be useful for 
signal and background discrimination. 

The discrimination distributions of multivariate analysis are demonstrated in Fig. 8. 
As it is shown, multivariate analysis methods (Neural Network and Boosted Decision Tree) 
using the input kinematic observables indeed work well to separate the signal from back- 
ground events. And the efficiencies of these two TMVA methods are also listed in the 
Tables 5 and 6, as a comparison to that of the simple cut method used above. As shown 
in the Tables 5 and 6, both of the neural network and boosted decision tree methods can 
improve the significance by a factor 100% or more. For instance, for the charged Higgs with 
mass 1 TeV, both TMVA methods can improve the significance by gaining more signals at 
the price for a little bit more background events. This improvement is also demonstrated 
in the mass window of the charged Higgs boson, as demonstrated in Fig. 9. 



I MLP Discriminant Distribution ] [ BDT Discriminant Distribution ] 




Figure 8: The discriminant distributions of neural network analysis and boosted decision tree 
analysis are demonstrated, where both signal and background is normalized to one. 

Sensitivity of LHC to heavy charged Higgs boson in the mu±-a{pp — > H^t) plane are 
investigated by scanning the charged Higgs boson mass rrijj± from 800 GeV to 1.5 TeV with 
a step 100 GeV. For each Higgs boson mass, the cone parameter R in the JH-tagger method 
with the best selection efficiency for signal is taken. The sensitivity reachable at LHC with 
14 TeV collision energy with two b taggings is provided in Table 7 and is displayed in Fig. 
(10). Roughly speaking, the TMVA can improve the sensitivity by a factor 100%. In Fig. 
(10), we show bounds from both 100 fb _1 and 500 fb _1 integrated luminosity. Two sample 
points of MSSM (marked with a solid star and a solid circle, respectively) are also shown. 
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Figure 9: The reconstructed mass bump of charged Higgs boson in three analysis methods are 
demonstrated, where both signal and background have been normalized to one before analysis. 





signal 
m H ± = 1.0 TeV 


tt + jets 


QCD 

nj >A 


After simple cuts 


2.5% 


1.5 x 1(T 5 


1.5 x 10" 8 


After NN cut (NN > 0.6) 


5.5% 


2.0 x 1(T 5 


3.0 x 1(T 8 


After BDT cut (BDT > 0.5) 


5.7% 


2.1 x 10" 5 


3.2 x 10~ 8 



Table 5: Comparison among the efficiencies of the neural network, boosted decision tree methods 
and hybrid-R reconstruction method by step-by-step cuts for M H ± = 1 TeV. 

It is found that the signal and background separation can be improved when the mass 
of the charged Higgs boson increases. From Fig. 10, it is reliable to conclude that the 
whole mass range of heavy charged Higgs boson in our study can be feasible at the LHC. 

5. Conclusions 

In this paper, we have studied the full hadronic decay mode of the heavy charged Higgs 
production process pp — > t2H~ — > t2^i^3- The main task to reconstruct the heavy charged 
Higgs boson in association with a non-boosted top comes from the large combinatorics and 
the huge QCD background events. By using the b taggings and top tagger methods to 
suppress the background from the SM, we propose a hybrid-R reconstruction method to 
fully reconstruct the signal. 

Our hybrid-R reconstruction method exploits the advantages of a larger and a smaller 
cone size of jet algorithms and can greatly reduce the combinatorics. As demonstrated 
in our analysis, this method is helpful to confirm the signature as well as to separate the 
signal from the background. 

A set of delicate kinematic observables are applied to the multi-variate analysis. Com- 
pared with the simple cut method, an improvement around 100% in the S/\fB can be 
achieved. The sensitivity of the charged Higgs boson mass m H ± from 800 GeV to 1.5 TeV 
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signal 
m H ± = 1.5 TeV 


tt + jets 


QCD 

»i>4 


After simple cuts 


1.8% 


1.9 x 10~ 6 


1.8 x 10~ 10 


After NN cut (NN > 0.7) 


4.0% 


2.0 x 10" 6 


3.6 x 10" 10 


After BDT cut (BDT > 0.5) 


4.5% 


2.1 x 10" 6 


2.4 x 10" 10 



Table 6: Comparison among the efficiencies of the neural network, boosted decision tree methods 
and hybrid-R reconstruction method by step-by-step cuts for M H ± =1.5 TeV. 



m H ± (TeV) 


0.8 


0.9 


1 


1.1 


1.2 


1.3 


1.4 


1.5 


a(fb) 


324 


192 


118 


75 


49 


33 


22 


15 


^= (with two b taggings & TMVA) 


19.4 


13.3 


9.5 


7.0 


5.1 


4.0 


3.1 


2.6 


lower bound on cr(fb) 


33 


29 


25 


21 


19 


16 


14 


11 


-4= (with two b taggings without TMVA) 


9.3 


6.5 


4.6 


3.4 


2.5 


1.9 


1.5 


1.3 


lower bound on cr(fb) 


70 


59 


51 


44 


39 


35 


29 


23 



Table 7 : The significance and sensitivity of LHC for heavy charged Higgs boson production are 
shown, where wc assume the total integrated luminosity as 100 fb _1 . For the sake of comparison, 
the results with and without TMVA are provided. We use ^= — 2.0 to obtain the exclusion lower 
bound on the cross section. 

is also explored and it is found that the LHC can probe the heavy charged Higgs boson 
mass up to 1.5 TeV. When the integrated luminosity is higher, even heavier mass region 
of the charged Higgs boson can be reachable. 

In this paper, we chose a reference THDM model as an example to perform a concrete 
study and the cross section of heavy charged Higgs boson is approximately proportional 
to Yl + Yft. To apply our results to other models, like the supersymmetry model, one can 
rescale the cross section and the significance in term of the model parameters. 

In the evaluation of cross section of signal, we have included the process gb — > btH~ 
which is a part of QCD next-to-leading (NLO) order correction of process gb — > tH~and 
we only assume the K- factor as 1.5 here. A lot of works on high order QCD corrections and 
SUSY corrections for process bg —> tH~ have been done in the past, including QCD NLO 
calculations [66, 67], SUSY-QCD NLO calculations [67], Yukawa and SUSY electroweak 
corrections [68] as well as high-order soft gluon corrections [69, 70, 71]. These calculations 
show the reduction of scale dependence and a significant enhancement to lowest order result 
when high order corrections are considered. For fj, = m#± , tan/3 = 30 and m^± = 800GeV 
(1 TeV), the K- factor in QCD next-to-next-to-next-to leading order soft gluon threshold 
correction at next-to-leading logarithmic accuracy is 1.75 (1.81) [70, 71]. Therefore when 
our results are applied to other models, the significance can be better. 

Obviously, our hybrid-R reconstruction method can be generalized to cases where some 
jets are decay products of one highly boosted massive particle while some of others come 
from a massive particle near the threshold. Such cases can exist in some new physics 
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Sensitivity to a heavy charged Higgs boson at LHC 
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Figure 10: The sensitivity curves of LHC to the heavy charged Higgs bosons with two b tagging 
are shown from both 100 fb _1 and 500 fb _1 datasets. The solid star and solid circle markers as two 
example points of MSSM are also shown. 

model where the production of a new heavy particles X in association with a SM massive 
particle Y, where X decays to boosted massive particles of the SM, such as W bosons, 
Z bosons, H bosons or top quarks. For examples, the hybrid-R reconstruction method 
can be applied to the processes pp — > W H — > tbH and pp — > t W — >• WbW (where W 
and t represent a new gauge boson and a new up-type quark, respectively). Additionally, 
our analysis can be extended to the semi-leptonic and di-leptoinc modes [72] and to the 
intermediate heavy charged Higgs boson region (say from 300 GeV to 800 GeV) with the 
help of HEPToptagger. 
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